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Abstract 

In this paper we prove that the lack of uniqueness for solutions of 
the tree dyadic model of turbulence is overcome with the introduction 
of a suitable noise. The uniqueness is a weak probabilistic uniqueness 
for all Z 2 -initial conditions and is proven using a technique relying on the 
properties of the q-matrix associated to a continuous time Markov chain. 



1 Introduction 

The deterministic dyadic model on a tree 



X Q (t) = 0, 



(1) 



was introduced as a wavelet description of Euler equations in [15] and stud- 
ied in [4J as a model for energy cascade in turbulence. It can be seen as a 
generalization with more structure of the so called dyadic model of turbulence, 
studied in [fj- As we show in section [5] this deterministic model ([I]) does not 
have uniqueness in I 2 . The aim of this paper is to prove that we can restore 
uniqueness with the introduction of a suitable random noise: 

dXj = ( Cj X? - c k XjX k )dt + aCjXj o dW 3 - a ^ c kX k ° dWk, (2) 
keOj keOj 

with (Wj)j e j a sequence of independent Brownian motions. Let's also assume 
deterministic initial conditions for $Z§: X (0) = x = (xj)- eJ £ I 2 . The main 
result of this paper is the weak uniqueness of solution for ([2]) , proven in theo- 
rem Ol 

This paper can be seen as a generalization to the dyadic tree model of the 
results proven for the classic dyadic model in [5] , but the proof of uniqueness 
given here relies of a new, different approach (see also jS]) based on a general 
abstract property instead of a trick (see Section [6]). The q— matrix we rely on 
is closely related to an infinitesimal generator, so the technique is valid for a 
larger class of models. 

The set J is a countable set and its elements are called nodes. We assume 
for the nodes a tree-like structure, where given j £ J, j is the (unique) father 
of the node j, and Oj C J is the finite set of offsprings of j. In J we identify a 
special node, called root and denoted by 0. It has no parent inside J, but with 
slight notation abuse we will nevertheless use the symbol when needed. 
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We see the nodes as eddies of different sizes, that split and transfer their 
kinetic energies to smaller eddies along the tree. To formalize this idea we 
consider the eddies as belonging to discrete levels, called generations, defined as 
follows. For all j G J we define the generation number \j\ £ N such that |0| = 
and |fc| = \j\ + 1 for all k E Oj. 

To every eddy j G J we associate an intensity Xj(t) at time t, such that 
Xj(t) is the kinetic energy of the eddy j at time t. The relations among intensi- 
ties are those given in ([1]) for the deterministic model and © for the perturbed 
stochastic model. The coefficients Cj are positive real numbers that represent 
the speed of the energy flow on the tree. 

The idea of a stochastic perturbation of a deterministic model is well estab- 
lished in the literature, see [5] for the classical dyadic model, but also [IT], [9] 
for different models. This stochastic dyadic model falls in the family of shell 
stochastic models. Deterministic shell models have been studied extensively 
in [12J while stochastic versions have been investigated for example in |10| and 



When dealing with uniqueness of solutions in stochastic shell models, the 
inviscid case we study is more difficult than the viscous one, since the more 
regular the space is, the simpler the proof and the operator associated to the 
viscous system regularizes, see for example [3] about GOY models, where the 
results are proven only in the viscous case. 

In @ the parameter a ^ is inserted just to stress the open problem of the 
zero noise limit, for er — > 0. This has provided an interesting selection result 
for simple examples of linear transport equations (see [2] ) , but it is nontrivial 
in our nonlinear setting, due to the singularity that arises with the Girsanov 
transform, for example in ((5]). 

It is worth noting that the form of the noise is unexpected: one could think 
that the stochastic part would mirror the deterministic one, which is not the 
case here, since there is a j-indexed Brownian component where we'd expect a 
j one, and there is a fc-indexed one instead of a j one. 

One could argue that this is not the only possible choice for the random per- 
turbation. On one hand we chose a multiplicative noise, instead of an additive 
one, but this is due to technical reasons (see [H]). On the other hand, there 
are other possible choices, for example the Brownian motion could depend on 
the father and not on the node itself, so that brothers would share the same 
Brownian motion. But the choice we made is dictated by the fact that we'd like 
to have a formal conservation of the energy, as we have in the deterministic case 
(see [4]). If we use ltd formula to calculate 



= (cjX^Xj - c k X)X k )dt + o-CjXjXj o dWj - a ^ c k X 3 X k o dW k , 



[16j. 




we can sum formally on the first n + 1 generations, taking Xq(€) = 0: 




b'l=o 



\j\—n k£Oj 



Y CkXjXkiXjdt + aodWk), 



\ j\—n k£Oj 
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since the series is telescoping in both the drift and the diffusion parts indepen- 
dently. That means we have P-a.s. the formal conservation of energy, if we 
define the energy as 

E n (t) = V X 2 (t) £(t) = VX 2 (t) = lim £ n (t). 
\j\<n je.J 

2 Non-uniqueness in the deterministic case 

In [7] it has been proven that there exists examples of non uniqueness of I 2 
solutions for the dyadic model if we consider solutions of the form Y n (t) — 

— — , called self-similar solutions, with (a n ) n € I 2 . Thanks to the lifting result 

t — to 

(Proposition 4.2 in [3]) that is enough to obtain two different solutions of the 
dyadic tree model, with the same initial conditions. 

Following the same idea of self-similar solutions, introduced in [7\ and [5J for 
the classic dyadic model and in section 5.1 in [3] for the tree dyadic model, we 
can construct a direct counterexample to uniqueness of solutions. In order to 
do this we need an existence result stronger than the one proven in [4]. 

Theorem 2.1. For every x € I 2 there exists at least one finite energy solution 
of with initial conditions X(0) = x and such that 

J2xf(t)<J2x 2 ( S ) V0<s<t. 
jeJ jeJ 

The proof of this theorem is classical, via Galerkin approximations, and 
follows that of theorem 3.3 in 0]. 

Now we recall the time reversing technique. We may consider the system ([T]) 
for t < 0: given a solution X(t) of this system for t > 0, we can define X(t) = 
— X{— t), which is a solution for t < 0, since 

j t X 3 {t) - jXj{-t) - Cj X 2 (-t) - £ c k X^t)X k {-t) 

keOj 

= c 3 X 2 (t)- ]T CkX^Xkit). 

We can now consider the self similar solutions for the tree dyadic model, 

as introduced in [4], Xj(t) = — - — , defined for t > to, with to < and with 

t — t 

(aj)j£j £ I 2 such that 

(a Q = 

S dj + Cja 2 - = ^2 Cka 3 a k , Vj £ J. 
I keo 3 

We time-reverse them and we define 

Xj(t) = -Xj(-t) Vj £ J t < -t , 
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which, as we pointed out earlier, is a solution of ([T]) in (—00, —to), with —to > 0. 
Since 

lim \Xj(t)\ = and lim \XAt)\ = +00, Vj £ J 

we have 

lim \XAt)\ = and lim \XAt)\ = +00, Mj e J. 

Thanks to theorem 12.11 there is a solution X, with initial conditions x = X(0), 
and this solution is a finite energy one, so, in particular, doesn't blow up in — <o- 
Yet it has the same initial conditions of X, so we can conclude that there is no 
uniqueness of solutions in the deterministic case. 

3 ltd formulation 

Let's write the infinite dimensional system ([2]) in Ito formulation: 

dXj = (cjXS - c k X 3 X k )dt + acjXjdWj 
keOj 

2 

CkX k dW k - °— (c 2 + cl)Xjdt. (3) 

We will use this formulation since it's easier to handle the calculations, while 
all results can also be stated in the Stratonovich formulation. 

So let's now introduce the definition of weak solution. A filtered proba- 
bility space (O, Tt, P) is a probability space (O, J-oo,P) together with a right- 
continuous filtration (J-"t)t>o such that Too is the cr-algebra generated by (J 00 ^t* 

Definition 1. Given x e I 2 , a weak solution of in I 2 is a filtered proba- 
bility space (O, J- t ,P), a J-indexed sequence of independent Brownian motions 
(Wj)j£j on (Q, J-"t,P) and an Z 2 -valued process (Xj)j e j on (O, J-"t,P) with con- 
tinuous adapted components X,- such that 



- (\ Cj X 2 {s)- CkX^X^da 

+ f aCjXj^dWjis) - V / ac k X k (s)dW k (s) 



y f ( c , 2 + E c *)^( s ) ds > ( 4 ) 



keOj 

for every j g J, with Co = and -X"o(i) = 0. We will denote this solution by 

(n,r t ,p,w,x), 

or simply by X. 
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Definition 2. A weak solution is an energy controlled solution if it is a solution 
as in Definition Q] and it satisfies 

<E*?) = i, 

jeJ jeJ 

for all t > 0. 

Theorem 3.1. There exists an energy controlled solution to §3§ in L°°(fl x 
{0,T},1 2 ) for(x 3 )el 2 . 

We will give a proof of this Theorem at the end of Section [7] It is a weak 
existence result and uses the Girsanov transform. 

We'll prove in the following result that a process satisfying Q satisfies ^ 
too. 

Proposition 3.2. If X is a weak solution, for every j 6 J the process (Xj(t)) t >o 
is a continuous semimartingale, so the following equalities hold: 

f t ac j Xj(s)odW j (s)= f acjX^dW^s)-^- f c^X^ds 
o Jo z Jo 

[ a Y c k X k (s) o dW k (s) = Y I <rc k X k ( S )dW k (s) + ^ f ^ c 2 k X j (s)ds, 

where the Stratonovich integrals are well defined. So X satisfies the Stratonovich 
formulation of the problem ^ . 

Proof. We know that 

ft ft 

nr. XJs)d.W,(s) 4- - ^ 



CTC " 

<TCjXj{8) o dW,{s) = / acjX-MdWiia) + -±[Xj, Wj]t, 



but from (f2]) we have that the only contribution to [Xj, Wj] is given by the 
—aCjXj o dWj term, so 

^[X 7 , Wi] t = ^-\-f ° c i X i ° dW^W^t = J* Xjds. 

Now if we consider the other integral, we have 

pt pt 

<rJ2 CkX k (s) o dW k (s) = Y / o-c k X k (s)dW k (s) + ^ -±[X k ,W k ] t . 

For each X k we get, with the same computations, that the only contribution to 
[X k , W k ]t comes from the term ac k Xj o dW kl so that we get 

^[X k , W k } t = ^[[ ac k X 3 o dW k , W k ] t = ^ / Xjds. □ 
z * Jo * Jo 
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4 Girsanov transform 

Let's consider © and rewrite it as 
dXj = CjXjiXjdt + adWj)- 

c k X k (X j( lt + adW k ) - G —{c) + 4) X jdt. (5) 

The idea is to isolate Xjdt + adWj and prove through Girsanov's theorem that 
they are Brownian motions with respect to a new measure P in (CI, Too), simul- 
taneously for every j G J. This way ([3]) becomes a system of linear SDEs under 
the new measure P. The infinite dimensional version of Girsanov's theorem can 
be found in [T7] and [T3]. 

Remark 1. We can obtain the same result under Stratonovich formulation. 

Let X be an energy controlled solution: its energy £(t) is bounded, so we 
can define the process 



M t = -i^ / Xj^dW^s) !(.) 
which is a martingale. Its quadratic variation is 

Because of the same boundedness of £ (t) stated above, by the Novikov criterion 
exp(M t — i[M] t ) is a (strictly) positive martingale. We now define P on (f2, Ft) 
as 



dP 
dP 



= exp(M t -hM] t ) 

Ft £ 



exp(-i]T / X^dWM ~ A [ ^2x§(s)ds) 



(7) 



for every t > 0. P and P are equivalent on each T t , because of the strict 
positivity of the exponential. 

We can now prove the following: 

Theorem 4.1. If (Q, J- t , P,W, X) is an energy controlled solution of the non- 
linear equation Q, then (17, J- t , P, B, X) satisfies the linear equation 

2 

dXj = acjXjdBjit) -<? c kX k dB k (t) - y [c] + ^ c t) x jdt, (8) 
keOj keOj 

where the processes 

B j (t)=W j (t)+ I -Xj( s )ds 
Jo o 

are a sequence of independent Brownian motions on (O, J-j,P), with P defined 
by 0. 



G 



Proof. Now let's define 



B j (t)=W j (t)+ [ lxj(s)ds. 



o 

Under P, {Bj (£))je./.te[o,T] is a sequence of independent Brownian motions. 
Since 

o~ f CjXj^dBjis) = a ( CjXj(s)dWj(s)+ [ c J X 1 2 (s)ds 
Jo Jo Jo 

a [ c k X k (s)dB k (s) = a [ c k X k (s)dW k (s) + [ c k Xj(s)X k (s)ds keOj. 
Jo Jo Jo 



Then (O can be rewritten in integral form as 



Xj(t) = x 3 +cr [ Cj Xj(s)dBj(s) - a V f c k X k {s)dB k {s) 

T A c i + £ c l) x ^ s ) ds ' ( 9 ) 
Jo keOj 

which is a linear stochastic equation. □ 
Remark 2. We can write our linear equation © also in Stratonovich form: 

dXj — acjXjo dBj(t) — ac k X k o dB k {t). 
keOj 

Remark 3. If we look at @ we can see that it is possible to drop the cr, 
considering it a part of the coefficients Cj . 

We can use Ito formula to calculate 

\dX] = XjdXi + \d[Xi\ t (10) 

= acjXjXjdBj — a S ' c k XjX k dB k 

keo-j 

~ y(s 2 + £ 4)x]dt+^(c 2 jX § + J2 clxl)dt 

keOj keOj 

= - y(s 2 + £ 4)*?* + dty + y(s 2 4 + E 4*2)*, 

keOj keOj 

Nj(t)=af CjXjXjdBj - a V / c k XjX k dB k . 

This equality will be useful in the following. 

We now present an existence result also for system ©. 

Proposition 4.2. There exists a solution of © m T°°(fi x [0,T],/ 2 ) wtt 
continuous components, with initial conditions x G I 2 . 



with 



(11) 
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Proof. Fix N > 1 and consider the finite dimensional stochastic linear system 

'dXf = aCjXfdBjit) ~0-J2 keO] c k X^dB k {t) 

~4(S 2 + T.ueo, 4)X?dt j G J, < |i| < N 

Xj*(t) = k e J, \k\ = N + 1 

xf(o)= Xj jeJ,o<\j\<N. 

This system has a unique global strong solution (X^)j e j. We can compute, 
using ([TU]) and the definition of iVj in (fTTjl , 



\j\<N \i\<N keOj 



- E t E ) a * 

y E 4(^f) 2 <0. 

\k\=N+l 



Hence 



(xf(t)) 2 < e ^ - a - s - vt > o. 



|i|<7V |j|<jv jeJ 



This implies that there exists a sequence iV m f oo such that (Xj Vm ) J - 6 j converges 
weakly to some (Xj)j e j in L 2 (£l x [0,T],/ 2 ) and also weakly star in L°°(£l x 
[0,T],; 2 ), so (Xj) je j is in L°°(f7 x [0,T],Z 2 ). 

Now for every JV e N, (Xj)j^j is inProg, the subspace of progressively 
measurable processes in L 2 (Q, x [0,T],l 2 ). But Prog is strongly closed, hence 
weakly closed, so (Xj)j e j £ Prog. 

We just have to prove that (Xj)j^j solves (JSJ. All the one dimensional 
stochastic integrals that appear in each equation in (0) are linear strongly con- 
tinuous operators Prog — > L 2 (f2), hence weakly continuous. Then we can pass 
to the weak limit in (fT!?)) . Moreover from the integral equations © we have 
that there is a modification of the solution which is continuous in all the com- 
ponents. □ 

5 Closed equation for Ep[Xj(t)] 



Proposition 5.1. For every energy controlled solution X of the nonlinear equa- 



tion @ 7 Ep[X 2 (t)] is finite for every j £ J and satisfies 



jEp[X 2 (t)] = -a\c) + £ 4)Ep[X 2 (t)} 

keOj 



a 2 c 2 E p [X 2 {t)] + a 2 E c 2 E p [X 2 (t)}. (13) 



keOj 



Proof. Let (Q, J-j, P, W, X) be an energy controlled solution of the nonlinear 
equation ([3]), with initial condition X £ I 2 and let P be the measure given by 
Theorem 14. II Denote by Ep the expectation with respect to P in (fi, .F t ). 
Notice that 

%[/ X/(t)dt]< oo VjGJ. (14) 
Jo 

For energy controlled solutions from the definition we have that P-a.s. 



^^(t)<max4(t)^X|(t)<(^ : 



a; 2 ) 2 



because of the behavior of the energy we showed. But on every J- tl P ~ P, so 

je.7 je.7 

and (fill) holds. 

From lfl"4"]) it follows that Mj(t) is a martingale for every j S J. Moreover 

£pE^(t)]<», 

since (i) is an energy controlled solution and the condition is invariant under 
the change of measure P <->■ P on JF t and, in particular, 

Ep[X 2 (t)} < oo Vj e J. 
Now let's write (fTU)) in integral form: 

X 2 (<) - x) = -a 2 /V 2 + J] c 2 fe )X 2 ( S )d S 

+ 2 / dN 3 (s) + a 2 f\c 2 X 2 ( S ) + £ 4X 2 k {s))d S . 



We can take the P expectation, 



Ep[X 2 {t)] - x 2 = -a 2 f\c 2 + £ c 2 k )E p [X 2 (s)}d S 

J ° k£0 3 
ft 



a 2 f c 2 Ep[X 2 ( S )}ds + a 2 £ / c 2 Ep[X 2 ( S )]d S , 



feeOj 

where the Nj term vanishes, since it's a P-martingale. Now we can derive and 
the proposition is established. □ 

It's worth stressing that Ep[X 2 (t)] satisfies a closed equation. Even more 
interesting is the fact that this is the forward equation of a continuous-time 
Markov chain, as we will see in the following section. 
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6 Associated Markov chain 



We want to show and use this characterization of the second moments equation 
as the forward equation of a Markov chain, taking advantage of some known 
results in the Markov chains theory. We follow the transition functions approach 
to continuous times Markov chains; we don't assume any knowledge of this 
theory, so we will provide the basic definitions and results we need. More results 
can be found in the literature, see for example [T]. 

Definition 3. A non-negative function fj.i(t) with j, I G J and t > is a 

transition function on J if fji(0) = 8ji, 

^/iK*)<i VjeJ, vt>o, 

lEJ 

and it satisfies the semigroup property (or Chapman-Kolmogorov equation) 
(i + s) = £ f jh (t)fhi (s) Vj, I G J, V*, s > 0. 

Definition 4. A q-matrix Q = {qji)j,i&j is a square matrix such that 
< q 3 i < +oo Vj ^ / G J, 
X! - = : Qj < +oo Vj G J. 

A q-matrix is called stable if all q^s are finite, and conservative if 

If Q is a q-matrix, a Q-function is a transition function fji(t) such that 

fji(o) = Q- 

The q-matrix shows a close resemblance to the infinitesimal generator of the 
transition function, but they differ, since the former doesn't determine a unique 
transition function, while the latter does. Still this approach can be seen as a 
generator approach to Markov chains in continuous times. 

Now let's see these objects in our framework: let's write (fT3|) in matrix form. 
Let Q be the infinite dimensional matrix which entries are defined as 

Qj,j = -°" 2 ( c i + c ?') q M = 0-2 c i q i' k = 1 {keo ] }0 2 c 2 k for k j= j, j 

keo-j 

Proposition 6.1. The infinite matrix Q defined above is the stable and conser- 
vative q-matrix. Moreover Q is symmetric. 

Proof. It's easy to check that Q is a stable and conservative q-matrix. First of 
all qj.j < for all j € J and qjj > for all j ^ I. Then 

'/./ = X! = + = a ~ r j + X! ff2c fe = 
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Moreover it is very easy to check that the matrix is symmetric: 



o z <A 1 = 3 jeO, 



*3 ' — \ 2 2 i - /0 ■ 7 2 2 f qt i 




□ 



Since Q is a g-matrix we can construct the process associated, as a jump and 
hold process on the space state, which in our case is the tree of the dyadic model. 
The process will wait in node j for an exponential time of parameter qj, and 
then will jump to j or k £ Oj with probabilities Qjj/lj and qj,k/lj respectively. 
This process is a continuous time Markov chain that has J as a state space and 
also has the same skeleton as the dyadic tree model, meaning that the transition 
probabilities are non-zero only if one of the nodes is the father of the other one. 

Given a g-matrix Q, it is naturally associated with two (systems of) differ- 
ential equations: 

VjlO = ^2vjh{t)QM (15) 

y'oitt) = ^IjhVhlit), 

called forward and backwards Kolmogorov equations, respectively. 

Lemma 6.2. Given a stable, symmetric and conservative q-matrix Q, then 
the unique nonnegative solution of the forward equations (|15p in L°°([0, oo), I 1 ), 
given a null initial condition y(0) — 0, is y(t) = 0. 

Proof. Let y be a generic solution, then 

jeJ 
2/i(0) = jeJ 

^2yj(t) < +oo. 

jeJ 

We can consider for every node yj — / + °° e~ t yj(t)dt 7 the Laplace transform 
in 1. From the last equation of the system above, we have y\ {jj < M, for some 
constant M > 0, so in particular we can consider k £ J such that y k > for 
all j G J. 

Now we want to show that y' k (t) is bounded: thanks to the symmetry and 
stability of Q we have 

\y'k(t)\ < | - q k yk{t)\ + | W(*)»fcl < % M + * M < +°°- 
We can integrate by parts 

r + co r+oo 

Vk = / ^ t y' k (t)dt= e~ ty ^yi{t)q lk dt = ^y iqik 

Jo Jo leJ leJ 

= -VkQk + ^2 yiqik - yk(~9k + X] 9kl) = 0, (17) 
l^k l^k 
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where the last equality follows from the conservativeness of Q, and we used the 
stability and symmetry. Now we have y^ — and so all yj = 0, hence yj(t) = 
for all j e J, for all t > 0. □ 

7 Uniqueness 

Now we can use the results of the previous section to prove the main results of 
this paper. 

Theorem 7.1. There is strong uniqueness for the linear system ((5|) in the class 
of energy controlled L°°(Q x [0,T},1 2 ) solutions. 

Proof. By linearity of © it is enough to prove that for null initial condi- 
tions there is no nontrivial solution. Since we have (|13[) . proposition 16.11 and 
lemma RT21 then Ep[Xj(t)] = for all j and t, hence X = a.s. □ 

Let's recall that we already proved an existence result for with proposi- 
tion! 



Theorem 7.2. There is uniqueness in law for the nonlinear system (J3J) in the 
class of energy controlled L 00 (Q x [0,T],/ 2 ) solutions. 

Proof. Assume that (fiW, jf\pM, W^,X^), i = 1, 2, are two solutions of © 
with the same initial conditions x G I 2 . Given n € N, t\, . . . , t n £ [0,T] and a 
measurable and bounded function / : (Z 2 ) n — »■ R, we want to prove that 

£ P(1) • ■ • (*»))] = E pW . . .,X^(t n ))}. (18) 

By theorem 14. II and the definition of P given in ((7]) we have that, for i = 1, 2, 

% (*«))] = 

£p w [exp{-M« + ±[M% MW] T }/(xW(t!), . . . ,*«(<„))], (19) 

where M^is defined as in ©. We have proven in proposition 14.21 and theo- 
rem 17.11 that the linear system @ has a unique strong solution. Thus it has 
uniqueness in law on C([0, T] , R) N by Yamada-Watanabe theorem, that is under 
the measures PW, the processes have the same laws. For a detailed proof 
of this theorem in infinite dimension see |17] , 

Now we can also include MW in the system and conclude that (jW, M^) 
under Pw have laws independent of i = 1, 2, hence, through (fT9"l) . we have (fT8")) . 

□ 

We can now conclude with the proof of Theorem 13.11 

Proof of Theorem\?n\ Let (0, T t ,P,B,X) be the solution of © in L°°(fi x 
[0, T] , Z 2 ) provided by theorem l7.ll We follow the same argument as in SectionHJ 
only from P to P. We construct P as a measure on (0, Ft) satisfying 



dP 
dP 



= exp(M T - - [M, M] T ) 
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where M t = ~ Y, je j £ x j( s ) dB j( s )- Under P the processes 

WAt) = B 3 (t) - / -X 7 (s)ds, 
Jo ° 

are a sequence of independent Brownian motions. Hence (17, JF t , P, W, X) is a 
solution of ([3]) and it is in L°°, since P and P are equivalent on Tt- □ 
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